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Abstract 

An exploratory study of chiral four-nucleon interactions in nuclear and neutron matter 
is performed. The leading-order terms arising from pion-exchange in combination with the 
chiral Air- vertex and the chiral NN37T- vertex are found to be very small. Their attractive 
contribution to the energy per particle stays below 0.6 MeV in magnitude for densities up 
to p = 0.4 fm -3 . We consider also the four-nucleon interaction induced by pion-exchange 
and twofold A-isobar excitation of nucleons. For most of the closed four-loop diagrams the 
occurring integrals over four Fermi spheres can either be solved analytically or reduced to 
easily manageable one- or two-parameter integrals. After summing the individually large 
contributions from 3-ring, 2-ring and 1-ring diagrams of alternating signs, one obtains at 
nuclear matter saturation density po = 0.16 fm -3 a moderate contribution of 2.35 MeV to 
the energy per particle. The curve E(p) rises rapidly with density, approximately with 
the third power of p. In pure neutron matter the analogous chiral four-body interactions 
lead, at the same density p n , to a repulsive contribution that is about half as strong. 
The present calculation indicates that long-range multi-nucleon forces, in particular those 
provided by the strongly coupled 7r./VA-system with its small mass-gap of 293 MeV, can 
still play an appreciable role for the equation of state of nuclear and neutron matter. 
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v©; 1 Introduction and summary 



According to their modern description and construction in chiral effective field theory, nuclear 
forces are organized in hierarchical way [TJ [2]. For generic few-body observables the contri- 
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q ■ butions arising from two-nucleon interactions are larger than those from three-nucleon forces, 
and the latter are again more important than possible corrections due to four-body forces. 
By constructing the chiral nucleon-nucleon potential up to next-to-next-to-next-leading order 

.^1 (N 3 LO) one has reached in the effective field theory approach the quality of a high-precision 
NN-potential in reproducing empirical NN-phase shifts and deuteron properties [U El [3] . By 
further lowering the resolution scale to A ~ 400 MeV the chiral NN-potential can be evolved to 
a low-momentum NN-potential Vi ow -k which becomes nearly model-independent and exhibits 
desirable convergence properties in perturbative calculations of many-nucleon systems and infi- 
nite nuclear matter [H El El E]- In chiral effective field theory the three-nucleon interaction can 
be constructed systematically and consistently together with the NN-potential [IJ|2]- At leading 
order it consists of a zero-range contact-term, a mid-range l7r-exchange component and a long- 
range 27r-exchange component, where the parameters of the latter occur also in the subleading 
27r-exchange NN-potential [£>1E]- The calculation of the subleading chiral three-nucleon force, 
built up by many pion-loop diagrams etc., has been completed in ref.[8]. The first calculation 
of the neutron matter equation of state with inclusion of these subleading chiral three-neutron 
interactions (at the Hartree-Fock level) has been presented recently in ref. [9]. 
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The leading four-nucleon interaction in chiral effective field theory has been constructed 
by Epelbaum in refs. [TOj HI] using the method of unitary transformations. The latter method 
allows to project the dynamics of the interacting pion-nucleon system into the purely nucleonic 
subspace relevant for few-nucleon systems below the pion-production threshold. The effect of 
the long-range part of the chiral four-nucleon force on the 4 He nucleus has been estimated in 
ref.p2]. An additional binding energy of the order of a few tenth MeV has been found by 
computing expectations values with totally symmetric wave-functions. Phenomenological four- 
nucleon forces arising from pion-exchange and twofold A-isobar excitation of nucleons have 
been considered in ref. [13] for computations of four-nucleon scattering processes (n- 3 H, p- 3 He, 
n- 3 He, p- 3 H and d-d). The effect of the four-nucleon force on the studied observables was found 
to be much smaller than that of the three-nucleon force, thus confirming the traditional belief 
in the hierarchical order of many-nucleon forces. However, the inclusion of the A-isobar was 
not able to resolve some longstanding discrepancies with experimental data [T3] . 

The present paper investigates the effects of chiral four-nucleon interactions on the equa- 
tion of state of homogeneous nuclear and neutron matter. The objective is to quantify those 
particular in-medium processes where four nucleons in the Fermi sea interact via the exchange 
of pions. The special case of pure neutron matter is interesting because the presence of neutral 
pions only leads to reduced isospin weight factors. The study of chiral four-body contributions 
to nuclear matter has been initiated in section 5 of ref . [2] , where a particularly simple class of 
2-ring Hartree diagrams (involving the chiral 47i~-vertex and the chiral NN37r-vertex) has been 
evaluated. With a contribution to the energy per particle of less than 0.1 MeV for densities up 
to twice normal nuclear matter density, 2po = 0.32 fm -3 , these long-range four-body correla- 
tions could be considered as negligibly small. In the present work we complete the leading order 
calculation by evaluating the corresponding Fock diagrams with one single closed nucleon ring. 
In addition to that we consider the four-nucleon interaction arising from pion-exchange and 
twofold excitation of nucleons to A-isobars. By introducing the direct coupling of the pion to 
the intermediate A-isobar (AA7r-vertex) a further type of long-range four-nucleon interaction 
can be generated. Although motivated more phenomenologically, these A-induced long-range 
four-nucleon interactions occur just as well in the chiral effective field theory framework. It is 
merely a question of how one counts the A- nucleon mass splitting (293 MeV), whether they 
are classified as leading order or as next-to-leading order terms. The calculations should be 
trustworthy for Fermi momenta sufficiently below the breakdown scale ~ 500 MeV of nuclear 
chiral effective field theory. 

The present paper is organized as follows: In section 2 we evaluate the contributions to 
the energy per particle of isospin-symmetric nuclear matter and pure neutron matter from the 
leading order four-nucleon interaction (involving the chiral 47r-vertex and the chiral NN37T- 
vertex). Let us note that our calculation does not follow the method of unitary transforma- 
tions of ref. [llj. In this scheme iterated two-body and three-body forces give rise to induced 
four-body forces, whereas we restrict ourselves to four-particle irreducible diagrams related 
to "genuine" four-nucleon forces (called V^ ass _ n in ref. [UJ). A complete study of the chiral 
four-body contributions to neutron matter following the method of unitary transformations 
has been performed recently by the Darmstadt-Ohio group [9]. Sections 3 and 4 deal with the 
A-induced four-nucleon interactions mediated by triple pion-exchange. The different orderings 
of the pion-couplings in the A-excitation-deexcitation process are combined into a symmetrized 
27r or 3tt contact-vertex proportional to the inverse AN mass-splitting or its square. Using this 
convenient short form we present results separately for 3-ring, 2-ring and 1-ring diagrams. In 
all except one case the occurring integrals over the product of four Fermi spheres can be either 
solved analytically or reduced to easily manageable one- or two-parameter integrals. As a major 
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result we find that the repulsive contribution from the dominant 3-ring diagram gets reduced to 
about 1/3 of its size by the (Fock-type) 2-ring and 1-ring diagrams. Altogether, there remains 
at saturation density p = 0.16 fm~ 3 a contribution of E(p ) = 2.35 MeV to the energy per par- 
ticle of isospin-symmetric nuclear matter. This moderate value lies on a curve for E(p) which 
however rises rapidly with density, approximately as p 3 . In pure neutron matter the analogous 
chiral four-body interactions lead, at the same density p n , to a repulsive contribution that is 
about half as strong. 

The present calculation indicates that long-range multi-nucleon forces, in particular those 
provided by the strongly coupled 7riVA-system with its small mass-gap of 293 MeV, are still of 
considerable size. Eventually, when combined with the relatively large attractive contributions 
from the subleading chiral 3N forces found in ref.[9] for pure neutron matter, one can expect 
cancellations to a certain degree. It remains as a future task to demonstrate these expected 
cancellations at a quantitative level for pure neutron matter and also for isospin-symmetric 
nuclear matter. 



2 Leading order four-body terms related to the chiral 

47r-vertex 

The leading long-range four-nucleon interaction is constructed by connecting the four nucleon 
lines by exchanged pions which couple according to the chiral pion-pion interaction. Since 
the off-shell chiral 47r-vertex is involved in the process this contribution alone depends on the 
choice of the interpolating pion-field [11] and thus is not unique. It has to be supplemented 
by the four-nucleon interaction generated additionally by the chiral NN37r-vertex, where the 
three pions emitted from one nucleon are absorbed on each of the other three nucleons. The 
sum of both pieces is actually representation-independent and thus gives rise to unique and 
physically meaningful results. The isovector (Weinberg- Tomozawa) NN7T7r-vertex is propor- 
tional to the energies of the exchanged pions, which in the present application are differences 
of nucleon kinetic energies. According to this property the four-nucleon interaction generated 
by two (Weinberg- Tomozawa) NN7T7T- vertices is a relativistic l/M 2 -correction. This observed 
suppression is equivalent to the statement Kiass-rn = in ref.[TT]. Note that we do not fol- 
low here the method of unitary transformations [11], but consider only "irreducible" diagrams, 
which would be commonly classified as "genuine" four-nucleon interactions. These pieces are 
called KLs-n in ref.[D]. 

Fig. 1 shows the diagrams one obtains by closing the four nucleon lines to either two rings 
or one ring. The combinatoric factors of these diagrams are 1/8, 1/2, 1/4 and 1, respectively. 
Diagram with more rings are trivially zero, due to a vanishing spin-trace. For the 2-ring 
(Hartree) diagrams the integrals over the four Fermi spheres factorize and can be readily solved 
by using the master formula eq.(29) in the appendix. One finds the following contribution to 
the energy per particle of isospin-symmetric nuclear matter: 
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with the dimensionless variable u = kf/m^. The nucleon density p is related to the Fermi 
momentum kf in the usual way, p = 2k 3 /3n 2 . The occurring parameters are: = 1-3 (nucleon 
axial- vector coupling constant), f n = 92.4 MeV (pion decay constant) and = 135 MeV 
(neutral pion mass). Note that after adding both 2-ring diagrams in Fig. 1 the only remainder 
of the chiral 7T7r-interaction is a constant factor —3m 2 /f 2 . This feature has ultimately lead to 
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Figure 1: 2-ring and 1-ring diagrams related to the leading order (genuine) chiral four-nucleon 
interaction. The combinatorial factor of each diagram is specified, if it is unequal to 1. 



the expression with a complete square in eq.(l). The analogous chiral four-body contribution 
in pure neutron matter is also of interest. In this case only neutral pions are present, which 
leads to a reduced isospin factor. Doing the calculation one finds from the sum of both 2-ring 
diagrams in Fig. 1 the following contribution to the energy per particle of pure neutron matter: 
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with u = kn/m^. The neutron density p n is related to the neutron Fermi momentum k n by 
pn = k\/2>it 2 . In order to simplify the notation, we make from here on the agreement that in 
all formulas for E n (p n ) the dimensionless variable u has the meaning u = k n /m n , while in all 
formulas for E(p) it keeps its original meaning u = kf/m n . 

The full and dashed line in Fig. 2 show these chiral four-body contributions to the equation 
of state of nuclear and neutron matter in the density region 0.04 fm -3 < p,p n < 0.32 fm -3 . 
One observes very small repulsive or attractive contributions which do not exceed 0.1 MeV in 
magnitude. Note that at equal densities, p n = p, the isospin reduction factor —1/6 for pure 
neutron matter gets effectively compensated by the increase of the Fermi momentum k n (by a 
factor 2 1 / 3 ). 

Next, we turn to the evaluation of the 1-ring (Fock) diagrams shown in the lower part of 
Fig. 1. In order to keep the symmetry of the left diagram (with a 4ir- vertex), it is advantageous 
to add four versions of the right diagram (with a NN37T- vertex) , each then weighted by a factor 
1/4. Putting together all pieces emerging from spin- and isospin-traces, one finds the following 
contribution to the energy per particle of isospin-symmetric nuclear matter: 
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Figure 2: Chiral four-body contributions to the energy per particle of nuclear and neutron 
matter arising from leading order 2-ring diagrams. 



with momentum transfers <fi = Pi — Pa, q 2 = P2 ~ Pi and = P3 — p 2 - 111 the case of pure 
neutron matter, where the chiral four-body interaction is mediated by neutral pions only, one 
gets from the sum of the 1-ring diagrams: 
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x 1 {qi qi-q 3 - 2q 2 ■ q 3 (qi + q 2 ) ■ <fi } . (4) 
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Note that we have exploited the symmetry under the interchange q\ -H- q 3 in order to simplify the 
expressions in the curly brackets of eqs.(3,4). The occurrence of a product of four different pion- 
propagators prohibits any substantial analytical reduction of these 12-dimensional integrals over 
the product of four Fermi spheres. Note also that there is a marked difference between E(p) 
in eq.(3) and E n (p n ) in eq.(4). After adding both 1-ring diagrams the result for pure neutron 
matter becomes proportional to m\, while there remains a finite term — 8[<fi • (q 2 x qs)] 2 for 
isospin-symmetric nuclear matter even in the chiral limit m n = 0. For comparison, the sum of 
the 2-ring (Hartree) diagrams has generated the m\ factor in both cases. We note as an aside 
that in the chiral limit, m % = 0, eq.(3) leads to the result E(p) ~ — 1.13 g\kj / (At: f^) 6 with a 
simple seventh power dependence on the Fermi momentum kf. 

The full and dashed line in Fig. 3 show the numerical results for the 1-ring (Fock) contribu- 
tions as a function of the density, p or p n . In the case of pure neutron matter one finds the usual 
pattern with the Fock contribution having opposite sign and reduced spin-weight in compari- 
son to the Hartree contribution (shown in Fig. 2). In contrast to this the Fock contribution for 
isospin-symmetric nuclear matter is about a factor of five larger (and of opposite sign) than the 
Hartree contribution. This reversed ordering comes from the fact that the Hartree contribution 
eq.(l) is "accidentally" suppressed by a factor (m n /kf) 2 , while the Fock contribution remains 
finite in the chiral limit m n = 0. The Hartree terms shown in Fig. 2 are therefore not represen- 
tative for the size of chiral four-body contributions. Irrespective of that exceptional behavior, 
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Figure 3: Chiral four-body contributions to the energy per particle of nuclear and neutron 
matter arising from leading order 1-ring diagrams. 

the attractive chiral four-body contributions represented by the diagrams in Fig. 1 stay in total 
below 0.6 MeV for densities up to p = 0.4 fm -3 . For this reason they can be considered as 
negligibly small. 

Let us remark that in the case of pure neutron matter the method of unitary transformations 
[TT] provides just one additional leading order chiral four-body interaction (called V a in ref.[9]). 
It has been computed in ref. with the result that its contribution to the energy per particle 
E n (p n ) reaches about —0.25 MeV at a density of p n = 0.2 fm -3 . 



3 Single delta-isobar excitation on two nucleons 

In this section we consider a more relevant type of long-range four-nucleon interaction [13] . It 
arises from triple pion-exchange and the excitation of two nucleons into virtual A-isobars. The 
corresponding four-nucleon interaction is depicted by the left diagram in Fig. 4. For an efficient 
treatment it is advantageous to combine the direct and crossed ttN — > A — > tcN transition 
into a contact-vertex for the absorbed and emitted pion. In the present application the energy 
dependence of the A-propagator i(u — A) -1 can be neglected since the energies to carried by the 
virtual pions are differences of small nucleon kinetic energies. With this valid approximation 
the pertinent contact-vertex reads: 

(5) 

where q a is an ingoing and % is an outgoing pion-momentum. We have inserted in eq. (5) already 
the empirically well-satisfied coupling constant ratio g^N^I 9-kNN = 3/v2, and A = 293 MeV 
denotes the mass-splitting between the A-isobar and the nucleon. Note that the scale A = 
293 MeV is comparable to the Fermi momentum kjo = 263 MeV at nuclear matter saturation 
density po = 0.16 fm -3 . At p n = 2po the neutron Fermi momentum k n = 418 MeV exceeds 
the scale A = 293 MeV already considerably. For the in-medium interaction the momenta of 
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Figure 4: Left: Long-range four-nucleon interaction with twofold A-isobar excitation. Right: 
Corresponding 3- ring Hartree diagram with combinatorial factor 1/2. 

the exchanged pions vary over the range < \q a ,b\ < 2/c/, n , but the included high- momentum 
part receives very little weight in the Fermi sphere integrals (as demonstrated in the appendix). 
When working with a NN7T7T contact-vertex the task of enumerating all topologically distinct 
in- medium diagrams gets simplified a lot. We turn now to the analytical evaluation of the 
pertinent four-loop in-medium diagrams representing the energy density, pE(p) or p n E n (p n ). 

3.1 Three-ring diagram 

This simplest closed in-medium diagram resulting from the A-induced four-nucleon interaction 
is the 3-ring diagram shown in the right part of Fig. 4. The upper two loops produce each a 
factor of density, p or p n , and the remaining integral over two Fermi spheres is readily solved 
by using eq.(29). One finds the following contributions to the energy per particle of isospin- 
symmetric nuclear matter: 

-. . g\mlv? [16m 6 in 4 20m 2 70m 3 n / 2 5\ W1 , 2 A fc . 

E{p) = (fkM^T - 12u + — + — arctan2M - [ 12u + 3) ln(1 + Au ) \ ' (6) 

and pure neutron matter: 

where the latter differs from the former only by an isospin reduction factor 1/12 (at equal Fermi 
momenta, k n = kf). 

Fig. 5 shows the four-body contributions to the energy per particle of nuclear and neutron 
matter as they arise from the 3-ring diagram with twofold A-excitation. For isospin-symmetric 
nuclear matter (see full line) one obtains a sizeable repulsion which reaches a value of 74 MeV at 
twice normal nuclear matter density 2p = 0.32 fm~ 3 . At saturation density po = 0.16 fm -3 the 
corresponding value E(p ) = 7.42 MeV is ten times smaller, but this amounts still to almost half 
of the empirical nuclear matter binding energy Eq ~ —16 MeV. The dashed curve in Fig. 5 for 
pure neutron matter lies about 20% below the full curve for isospin-symmetric nuclear matter. 
This feature demonstrates that the isospin reduction factor 1/12 in eq.(7) gets to a large extent 
compensated by the increase of the neutron Fermi momentum k n . 
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Figure 5: A-induced chiral four-body contributions to the energy per particle of nuclear and 
neutron matter: 3-ring diagram. 

With the analytical formulas eqs.(6,7) and the corresponding numerical results displayed 
in Fig. 5 the size of long-range A-induced four-nucleon contributions in nuclear and neutron 
matter is set. From here on the interesting question is, how much of a reduction will be provided 
by the Fock-type diagrams with fewer rings. 

3.2 Two-ring diagrams 

The 2-ring diagrams resulting from the chiral four-nucleon interaction with twofold A-isobar 
excitation are shown in Fig. 6. The last two diagrams (marked by 0) vanish identically. The 
upper diagram of this subset involves a central pion which carries zero-momentum and therefore 
does not couple according to the contact- vertex eq.(5). In order not to vanish in the first step, 
the spin-traces in the lower diagram select from the contact- vertex eq.(5) the spin-dependent 
part. Furthermore, by taking into account momentum conservation at each vertex one obtains 
a scalar triple product of three linearly dependent vectors (<fi + q2 + <?3 = 0), which vanishes 



Obviously, the upper two (topologically distinct) diagrams in Fig. 6 give rise to equal contri- 
butions. The left loop produces just a factor of density and the Fermi sphere integrals over the 
pion-propagators and momentum-dependent interactions on the right can be factorized with 
the help of tensors. Performing this strategy, as outlined in eq.(30) of the appendix, one ends 
up with the following expressions for the contributions to the energy particle of nuclear and 
neutron matter: 
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Figure 6: 2-ring diagrams resulting from the chiral four-nucleon interaction with twofold A- 
isobar excitation. 

Here, we have introduced the auxiliary functions: 
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defined by the tensorial Fermi sphere integral in eq.(30) of the appendix. The other auxiliary 
functions Hg T (x) occurring in eqs.(8,9) are closely related to Gs,t(x): 
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The lower left 2-ring diagram in Fig. 6 (with a combinatorial factor 1/2) is more difficult to 
evaluate. In order to facilitate an analytical treatment one chooses the momentum transfers 
carried by the exchanged pions as the basic variables and works with integrals over shifted 
Fermi spheres. In the end one arrives at a double-integral representation for the contributions 
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to the energy particle of nuclear and neutron matter: 
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By inspecting the first lines in eqs. (14,15) one realizes that the master integral eq.(29) has 
been employed. In addition to that, the following reduction formula has helped to eliminate a 
complicated angular integral: 



where the argument of F(. . .) on the left hand side gives the distance to the boundary of a shifted 
Fermi sphere, with p being the displacement of the center. We have checked numerically the 
analytical expressions in eqs. (14,15) against the original (untreated) four Fermi sphere integrals 
and found good agreement. The advantage of the double-integral representations for E(p) and 
E n (p n ) in eqs. (14, 15) is that these can be evaluated with much higher numerical precision. 

The full and dashed line in Fig. 7 show the A-induced chiral four-body contributions to the 
energy per particle of nuclear and neutron matter as they arise from the 2-ring diagrams. The 
first part given by eqs. (8,9) is typically a factor 3 to 4 larger than the second part given by 
eqs. (14, 15) and both pieces come with the same sign. By comparison with Fig. 5 one observes 
that the individually large contributions from the 3-ring and 2-ring diagrams nearly cancel 
each other, leaving a small repulsive (attractive) remainder for nuclear (neutron) matter. In 
magnitude this remainder is only about 1/7 of the starting values. As a consequence of this 
feature the final result for the A-induced chiral four-body contributions will be determined in 
a crucial way by the 1-ring diagrams. 

3.3 One-ring diagrams 

The 1-ring diagrams resulting from the chiral four-nucleon interaction with twofold A-isobar 
excitation are shown in Fig. 8. For the first diagram with a U-shaped pion-line three of the 
four occurring Fermi sphere integrals over the pion-propagators are structurally equivalent and 
with inclusion of the momentum-dependent interactions they factorize with the help of tensors. 
In the end only one single radial integral remains and the corresponding contributions to the 
energy per particle of nuclear and neutron matter read: 
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Figure 7: A-induced chiral four-body contributions to the energy per particle of nuclear and 
neutron matter: 2-ring diagrams. 




Figure 8: 1-ring diagrams resulting from the chiral four-nucleon interaction with twofold A- 
isobar excitation. 
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with the auxiliary functions Gs,t{x) given in eqs.(10,ll). The next two (topologically distinct) 
1-ring diagrams in Fig. 8 have a pion-line running zig-zag. Since both give equal contributions 
their combinatorial factor 1/2 gets effectively compensated. Again, the Fermi sphere integrals 
over the pion-propagators associated to the short pion-lines factorize via tensors. For the 
remaining integral over two Fermi spheres \pi$\ < kf the angular integrations can be carried 
out and in this procedure the third pion-propagator introduces the logarithmic function: 
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Putting all pieces together and exploiting the symmetry under the interchange x •<->■ y, one 
finds the following double-integral representations for the contributions of the zig-zag 1-ring 
diagrams to the energy per particle of nuclear and neutron matter: 
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The last 1-ring diagram in Fig. 8 involves crossed pion lines. The assignment of momenta to 
the virtual pions is now such that the previous factorization procedure does not work anymore. 
From the four Fermi sphere integrals only the one associated to the nucleon line between the 
contact-vertices can be performed analytically and the contributions to the energy per particle 
of nuclear and neutron matter read: 
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E n (Pn) = 
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(23) 



with momentum transfers q\ = f>\ — P3, q\ = P2 — P3 and the vector fj — (pi + p 2 — p 3 )/m n . 
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Figure 9: A-induced chiral four-body contributions to the energy per particle of nuclear and 
neutron matter: 1-ring diagrams. 



The full and dashed line in Fig. 9 show the A-induced chiral four-body contributions to 
the energy per particle of nuclear and neutron matter as they arise from the 1-ring diagrams. 
The three pieces (with U-shaped, zig-zag and crossed pion-lines) provide approximately equal 
amounts. As expected the sign of E(p) and E n {p„) is opposite to the contributions from the 
2-ring diagrams shown in Fig. 7. Although the magnitude of the repulsive 1-ring contributions 
is substantially reduced (by a factor 3 to 4) in comparison to the previous terms, these are by 
no means small. Such strongly with density rising four-body correlations, which reach values 
of 14.4 MeV or 21.2 MeV at 2p = 0.32 fm -3 , affect the equation of state of nuclear or neutron 
matter appreciably. 



4 Double delta-isobar excitation on one nucleons 

In this section we consider a further type of A-induced long-range four-nucleon interaction 
Once the A-isobar is excited, it can also couple directly to the pion (through the A Air- vertex). 
Each of the three pions emitted from the central nucleon gets absorbed on a spectator nucleon. 
For an efficient treatment of this chiral four-nucleon interaction it is advantageous to condense 
the three-pion dynamics at the central nucleon into a symmetrized three-pion contact-vertex 
as symbolized in Fig. 10. Approximating the A-propagator by the inverse mass-splitting —i/A 
and summing over all six permutations of (a, b, c) the pertinent NN37T contact-vertex reads: 



9l 



- 1 - 75e abc q a ■ (q b x q c ) + q a -q b a- q c (I85 ab r c - 75 ac r b - 75 bc r a ) 



40/3A2 

+q a -q c (185 ac r fe - 75 ab r c - 75 bc T a ) +q b -q c a-q a {I85 bc r a - 75 ac T b - 75 ab T c ) \ , (24) 



where q a , q b and q c denote outgoing pion- momenta. In order to fix the A An- vertex we use 
the coupling constant ratio g n AA/ 9ttNN = 1/5 of the quark-model. The relation for isospin 
(transition) operators, T a Q b T^ = {hi e abc — 6 ab r c + ^5 ac r b — 5 bc r a )/3, with 06 the 4x4 A-isospin 
matrices has been employed together with an analogous relation for spin (transition) operators. 
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Figure 10: The filled ellipse symbolizes the symmetrized three-pion contact-vertex involving 
double A-excitation. The associated chiral four-nucleon interaction gives rise to one non- 
vanishing 2-ring diagram and one non-vanishing 1-ring diagram. 



The 2-ring and 1-ring diagrams resulting from the three-pion contact-vertex are shown in 
Fig. 10. The upper 2-ring diagram vanishes because the spin-trace (on the left) selects the 
term proportional to the scalar triple product from eq.(24) with linearly dependent vectors 
{Qa + Qb + Qc — 0). For the lower 2-ring diagram (with combinatoric factor 1/2) the integrals 
over the Fermi sphere factorize and can be solved by using the master formula eq.(29). The 
corresponding contributions to the energy per particle of nuclear and neutron matter read: 
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The last 1-ring diagram in Fig. 10 can be treated in the same way as the zig-zag diagram in 
Fig. 8. One finds the following double- integral representation for its contribution to the energy 
per particle of nuclear and neutron matter: 
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(27) 
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Figure 11: A-induced chiral four-body contributions to the energy per particle of nuclear and 
neutron matter: 2-ring diagram. 
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with the logarithmic function L defined in eq.(19). Note that the peculiar numerical coeffi- 
cients entering the contact- vertex eq.(24), in particular the quark- model ratio 1/5, have all 
disappeared in the final expressions for E(p) and E n (p n ) in eqs. (25-28). 

Numerical results for the energy per particle of nuclear and neutron matter as they arise 
from the double A-induced long-range four-nucleon interaction are shown in Figs. 11,12. One 
observes attractive contributions from the 2-ring diagram which, in the next step, get nearly 
cancelled by repulsive contributions from the 1-ring diagram. After this balance the differ- 
ences between nuclear and neutron matter which are individually visible in Figs. 11,12 have 
essentially disappeared. The basic reason for the small net result in comparison to section 
3 is the absence of a 3-ring diagram (with its large spin-isospin weight factor) and the lower 
number of subsequent diagrams. The pattern of alternating and largely compensating 3-ring, 
2-ring and 1-ring contributions shown in Figs. 5,7,9,11,12 is specific for the three-pion exchange 
four-nucleon interaction with two-fold A-excitation considered here. When using fully anti- 
symmetrized many-nucleon forces [H [9] this merely technical decomposition into different ring 
diagrams does not come into play at all. 

5 Combined results and discussion 



After having evaluated various classes of long-range four-nucleon correlations in nuclear and 
neutron matter in the previous sections let us now draw a balance. The summed result of 
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Figure 12: A-induced chiral four-body contributions to the energy per particle of nuclear and 
neutron matter: 1-ring diagram. 




Figure 13: Chiral four-body contributions to the energy per particle of nuclear and neutron 
matter summed together. 



16 




density p or p n [fm ] 

Figure 14: Three-body contributions to the energy per particle of nuclear and neutron matter 
arising from two-pion-exchange with A-isobar excitation. 

all these contributions is shown in Fig. 13. For isospin-symmetric nuclear matter one finds at 
saturation density po — 0.16 fm -3 a moderate repulsive contribution to the energy per particle of 
2.35 MeV. However, the corresponding curve for E(p) rises strongly with the density and reaches 
the value 23.3 MeV at 2p = 0.32 fm~ 3 . One should note that in this density range, < p < 2p , 
a realistic nuclear matter equation of state resides in the binding regime with a negative energy 
per particle. A mere addition of long-range four-nucleon correlations as given by the full line in 
Fig. 13 would produce enormous distortions. In pure neutron matter the long-range four-body 
correlations come out about half as large (see dashed line in Fig. 13). To put it into proportion, 
a value of E n (2po) = 11.6 MeV amounts to about 1/4 of what sophisticated neutron matter 
calculations [15] give at that density. It is also interesting to compare with the three-body 
correlations arising from 27r-exchange with (single) A-excitation in the same framework. The 
corresponding analytical formulas for E(p) and E n (p n ) are written in eqs. (5,6,43,44) of ref. [16J 
and the associated numerical results are reproduced here in Fig. 14. By comparison with Fig. 13 
one observes that for isospin-symmetric nuclear matter these long-range three-body correlations 
are at 2po about a factor 1.5 larger than the four-body correlations, whereas in pure neutron 
matter the relative factor between both is about 2. Note also that three-body and four-body 
correlations grow approximately quadratically and cubically with the density, respectively. The 
hierarchical order of multi-nucleon forces is therefore manifest for the long-range components 
considered in this work. Obviously, it gets more pronounced the lower the density (p or p n ) is. 

In a recent work by the Darmstadt-Ohio group [9] a calculation of the neutron matter 
equation of state E n (p n ) with inclusion of the subleading chiral three-nucleon forces has been 
performed for the first time. These authors find relatively large attractive contributions, which 
reach about —4 MeV at the highest considered density p n = 0.2 fm~ 3 (see Fig. 2 in ref. [9]). This 
value is to be compared with 2.5 MeV repulsion due to chiral four-body correlations studied in 
the present work. It is therefore conceivable that substantial cancellations will occur among 
these novel higher order correlations. It remains as a future task to demonstrate the expected 
cancellations at a quantitative level for pure neutron matter and also for isospin-symmetric 
nuclear matter. 
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Figure 15: Weight functions for integrals over two and four Fermi spheres, respectively. The 
area under both curves is equal to 1. 

Appendix: Reduction of integrals over Fermi spheres 

In this appendix we collect several formulas which are useful for reducing integrals over (multi- 
ple) Fermi spheres. If the integrand depends only on the magnitude of the momentum transfer, 
the following reduction formula holds: 

|pi,2|<fc/ 

The weight function s 2 (l — s) 2 (2 + s) is slightly asymmetric about the midpoint s = 1/2 and it 
reaches its maximum at s = (y/A.2 — l)/2 = 0.5247. These features are exhibited by the dashed 
line in Fig. 15. In order to give a geometrical interpretation of this weight function, we note that 
the overlap volume of two unit-spheres with their centers displaced by 2s is 2n(l — s) 2 (2 + s)/3. 

Pions with their momentum-dependent interactions lead often to the following tensorial 
integral over a Fermi sphere of radius kf. 

f d 3 p 1 (p-pi)i(p-pi)j ml r 

J Wf ™* + (p-fl) 2 = ^f 5,J GS(X > U) + " %) Gt(X > ■ m 

Wi\<kf 

By contracting this equation with Sij and 3pipj — Sij the functions Gs(x,u) and Gt(x,u) 
can be calculated separately as simple scalar integrals over a Fermi sphere \p±\ < kf. The 
dimensionless variables on the right hand side are x = \p\/m n and u = kf/m w . The explicit 
analytical expressions for the functions Gs,t(x,u) involving arctangents and logarithms are 
given in eqs.(10,ll), suppressing in the notation the second argument u. 

Let us add a further remarkable reduction formula for integrals over four Fermi spheres. We 
assume a dependence of the integrand on the magnitude of the total momentum \pi+f>2+Pz+Pi\ 
only and derive the following reduction formula: 

r d l2 v 8k 12 r 1 

J j2^n\Pi+P2+n+P,\) = m ^2 Y J dsw 4 (s)F(4sk f ) , (31) 

Wj\<k f 



18 



with the (non-smooth) weight function: 

w 4 (s) = 4s(l-s) 6 (16s 4 + 96s 3 + 156s 2 + 56s-9) 

+6(1 - 2s) s(l - 2s) 6 (36 + 77s + 24s 2 - 12s 3 - 4s 4 ) . (32) 

This peculiar result for w^(s) has been obtained with Fourier transformation techniques and 
eq.(31) has been checked numerically for many examples of F(Askf). By employing the Fourier 
representation of the delta-function 5 3 (pi+p2+P3+P*4 — 4/cjs) the weight function w^(s) results 
with a multiplicative factor 7r/14175 from the integral s 2 / °° dq q 2 jo(4sg)[ji(g) / g] 4 , where jo,i(o) 
denote spherical Bessel functions. The weight function w^s) is normalized to J^dsw^s) = 
175/128 and a quick analysis shows that it reaches its maximum value of 3.63 at s = 0.335. 
The asymptotic behavior is: u>4(s) = 85s 2 for s — > 0, and u>4(s) = 1260(1 — s) 6 for s — > 1. 
These feature are exhibited by the full line in Fig. 15. The two weight functions put side by 
side in Fig. 15 demonstrate furthermore that multiple Fermi sphere integrals are dominated by 
total momenta of the order of the Fermi momentum kf. 
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